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Abstract. We explore the possibility of extending the well-known Berczin- 
Toeplitz quantization to reproducing kernel spaces of vector-valued functions. 
In physical terms, this can be interpreted as accommodating the internal de- 
grees of freedom of the quantized system. We analyze in particular the vector- 
valued analogues of the classical Segal-Bargmann space on the domain of all 
complex matrices and of all normal matrices, respectively, showing that for the 
former a semi-classical limit, in the traditional sense, does not exist, while for 
the latter only a certain subset of the quantized observables have a classical 
limit: in other words, in the semiclassical limit the internal degrees of freedom 
disappear, as they should. We expect that a similar situation prevails in much 
more general setups. 



1. Introduction 

Let n be a symplectic manifold, with symplectic form co, and H a subspace of 
L 2 (fl, dfi), for some measure /i, admitting a reproducing kernel K. For <j> £ C°°(f2), 
the Toeplitz operator with symbol 4> is the operator on TL defined by 

T+f = P(0/), fen, 

where P : L 2 (n,dfi) — > TL is the orthogonal projection. Using the reproducing 
kernel K, this can also be written as 



T*f(?)= / f(y)cf>(y)K(x,y)dti(y). 



n 

It is easily seen that is a bounded operator whenever <p is a bounded function, 
and HT^Ij-H^-H < ||0||oo, the supremum norm of <fi. 

Suppose now that both the measure /i and the reproducing kernel subspace TL 
are made to depend on an additional parameter h > (shortly to be interpreted as 
the Planck constant), in such a way that the associated Toeplitz operators Ti 1 ' on 
Hh satisfy, as h \ 0, 

(i-i) II^Hw^w,, - H\U 

and 

(1-2) WT^T^-T^Wn^n^Q, 

(A 1\ II 2ir \rp{h) rp{h)-i rp{h) n „ 
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(where {•, •} is the Poisson bracket with respect to to), and, more generally, 

oo 
3=0 

for some bilinear differential operators Cj : C°°(fi) x C°°(f2) — ► C°°(f2), with 
Co(</>, ■0) = 0V 1 an d Ci(0>V>) — Ci(tp,<j)) = ^r{0, "0}- Here the last asymptotic 
expansion means, more precisely, that 

N 

(1-4) \Tf ) Tf ) -Y J h j T^ ) \=0{h N+1 ) s S h\0, ViV = 0,1,2,.... 

i=o 

One then speaks of the Berezin-Toeplitz quantization. Indeed, it is well known that 
the recipe 

oo 

i=o 

then gives a star-product on f2, and i|1.2[l . (|1.3(l just amount to its correct semiclas- 
sical limit. 

Observe that if we introduce the normalized reproducing kernels — or coherent 
states — by 

K(;y) . , , , K(x,y) 



ky '- \\K(;y)\\n U - ky{x) ~ K(y,yy/*> 
and define the Berezin transform of an operator T on TL by 

T(y) := (Tk y ,k y ) n , 

then H1.4f) implies that as h — ► 0, 



J=0 



pointwise and even uniformly on f2. Often, one also has a stronger version of Ijl.lfl . 
namely 

(1.6) T^ l \x) -> <p(x) VxGO. 

The simplest instance of the above situation is f2 = R 2 ™ ~ C™, with the standard 
(Euclidean) symplectic structure, and 

(1-7) H h =L 2 hol ({l,dv h ) 

the Segal-Bargmann space of all holomorphic functions square-integrable with re- 
spect to the Gaussian measure dfih(z) := e~H z H / h (nh)~ n dz [dz being the Lebesgue 
measure on C"). The space Tih admits the reproducing kernel Kh(x, y) = e^ y ^ h . 
The Berezin transform turns out to be given just by the familiar heat operator 
semigroup, 

T£\x) = (irh)- n I e-^-yW 2 ' h <t>(y)d V , 
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so that, by the stationary phase (WJKB) expansion (see Section H3 below; for a 
discussion of the (WJKB) approximation, see for example |Me| . Chapter 7) 1 

(1.8) lf\x)=f2^ct>(x). 



=o J 



Similarly, 



T^T^\x) = (7rh)- 2n j j 0(y)V(2)e« a! ^+^>+^>-ll a: ll 2 -^ll 2 -ll z ll 2 )/ fe dy^, 
so by stationary phase again, 

(1.9) TfifKx) = dP+1 ^^ x ) h W\+\P\+\i\_ 

(Here the summation extends over all multiindices a, (3, 7, i.e. n-tuples of nonnega- 
tive integers a = {ol\, . . . , a n ), etc., and we are using the usual multiindex notations 
|a| =ai + •■ • + «„, a\ = a 1 \...a n \, d a = cH/ftc" 1 ...dx%».) Inserting JTBJ and 
(11.91) into l|1.5p. we get formulas for the cochains Cy. 

(1.10) C i(^)=Es fl °^ 5 ° i 

The resulting star-product coincides, essentially, with the familiar Moyal product. 

Other examples of Bcrezin-Toeplitz quantization include the unit disc D with 
the Poincare metric, bounded symmetric domains, strictly pseudoconvex domains 
with metrics having reasonable boundary behaviour, or, provided one allows not 
only holomorphic functions but also sections of line bundles as elements of Tth, all 
compact Kahler manifolds whose Kahler form is integral. In all these cases, the 
choice of the spaces (|1.7|) which works are the weighted Bergman spaces Tth = 
L^ ol (0, e - */^"), where n is the complex dimension of O and $ is a Kahler po- 
tential for u (so, for instance, for the unit disc Tth = L 2 wl (D, (1 — \z\ 2 Y 1 / h ^ 2 )). 
See |KS| . |BMS| or | AE| for the details and further discussion. 

In this note, we explore the possibility of extending the above formalism to 
reproducing kernel spaces of vector-valued functions. In physical terms, this can 
be interpreted as accommodating the internal degrees of freedom of the quantized 
system. 

In more concrete terms, this means that we again consider, for a given h > 0, 
a suitable measure /x/j on O, and a reproducing kernel subspace Tth C i^j, (f2, d[ih), 
where the subscript indicates that we are now dealing with vector-valued func- 
tions taking values in for some N > 1. The reproducing kernel Kh of Tth 
will thus now be a matrix-valued object, Kh : x f2 — > C NxN . We can again 
consider, for each <f> G C°°(f2), the associated Toeplitz operators, and investigate 
the existence of the asymptotic expansion 11.4J1 . In fact, we can now even allow 
matrix-valued symbols <f> £ C^ > JVX ]v(^)- By analogy with the scalar- valued case, 
one may again expect appropriate asymptotic expansions 
, — 00 

(1.11) T^ h \x) = J2 LrfWb*,, 

3=0 

^Throughout this paper, we are using the slightly nonstandard Laplacian A = £V d 2 /dzjdzj , 
which differs from the usual one by a factor of 4. 
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(1.12) T^\x)=Y,M j (4,,i;)(x)h j , 



with some differential and bidifferential operators Lj and Mj, respectively, whose 
comparison would yield l|1.5|) , thus suggesting that l|1.4|) is also likely to hold. 

Unfortunately, at this level of generality the results are negative: we show that for 
certain spaces Tih as above, which are quite natural generalizations of the Euclidean 
situation (|1.7|l to vector- valued functions, the semiclassical expansions (CP} i 
fail to hold. More specifically, it seems that there are asymptotic expansions for 

T^ h \x) and T^T^ l \x) in powers of h, but their coefficients do not depend only 
on the jets of </> and ip at x i but also at other points (i.e. are not local operators); 
besides, in addition to integer powers of h, half-integer powers seem to also enter 
the picture. Finally, l|l.lfl and l|1.6fl may also break down. 

However, it turns out that upon restricting to appropriate domains fi and ap- 
propriate classes of functions (f>, tp, the situation can be saved completely: namely, 
the following picture emerges. The admissible functions <fi can be identified with 
functions f(di; d,2, ■ ■ ■ , (In) on C x C^ -1 that are symmetric in the N — 1 variables 
d,2, . . . , c?at. Clearly, one can associate a function u# of this form to any function 
u : C — > C by the recipe u#(di; e?2, ■ ■ • , ^aO := u(di). We show that for any func- 
tions /, g of the above form, there exist uniquely determined functions u r on C, 
r = 0, 1, 2, . . . , such that 

(1.13) Tfri h) *f^h r T$. 

r=0 

Further, the u r are given by differential expressions involving / and g; and, finally, 
if / and g themselves are of the form v# and w# , respectively, for some functions 
v,w on C, then in fact u r — C r (v,w) with the bidifferential operators C r given 
by i|1.10fl for n = 1, i.e. £1 = C. This suggests the following interpretation: our 
quantum system has N internal degrees of freedom, which have no classical couter- 
parts, so that only a subset of the quantized observables have a classical limit. 
In the semiclassical limit the internal degrees of freedom disappear, as they should. 
We conjecture that similar quantizations can be carried out by our method in much 
more general setups. 

The paper is organized as follows. In Section |2 we introduce the spaces Tih 
of vector- valued functions on C 7VxAr , as well as their analogues for the subset of 
all normal matrices in C NxN . These spaces have previously appeared in AEG . 
In Section 13 we define a generalization of the Berezin transform — which will now 
also be a matrix-valued object — and establish its basic properties. Sections 01 
and [5] discuss the semiclassical asymptotic expansions <|l.lll) - (|1.12|) hi the above 
two settings of the domains of all matrices and all normal matrices, respectively. 
In Section [S] we present our first result in the positive direction, by establishing 
an asymptotic expansion — which is, however, of a highly non-local nature — for 

lf> and T^T^ in the case of the normal matrices. Finally, in Sections and [Sj 
we introduce our restricted class of observables <j) an( i establish the asymptotic 
expansion l|1.13f) . 
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2. The domains and the spaces 
Our first domain is f2 = C NxN , with the measures 

dn h {Z) =e -Wz)/h ^ h) -N* ^ 

where dZ denotes the Lebesgue measure on C NxN . The measures fXh are normal- 
ized to be of total mass one. The functions $fj(Z) := Z^ satisfy 

(2.1) / Z* 1 Z k d^i h {Z) = 5 ]k h k c k I 

Jn 

for some numbers Ck > 0; see AEG . Explicitly, ct are given by f |Gin| . formula 
(1.40), and ES|) 2 

(k + N + 1)1 



Ck 



that is, 



for k > N - 1, 
W(k + l)(k + 2) 

{k + N + iy. N\ 



N\(k+l)(k + 2) (k+l)(k + 2)(N -k-2)\ 



for k < N - 1, 



(2.2) c k 



k+i fc+i 



(fc + l)(fc + 2) 

In particular, Co = 1, ci = 2V, C2 = A^ 2 + 1, etc. 

It follows that if xi, . . . , \n is the standard basis of C^, then the functions 

(2.3) -4=' j = l,---,N, fc = 0,l } 2,..., 

are orthonormal in L 2 -,,, (f2, d/i^). Let be the subspace spanned by these func- 
tions. Then the function 

00 vk~y*k 

(2.4) ^ y ) = E^F- 

converges for all X,Y £ Q and is the reproducing kernel of Tin, in the sense that 

/ iv ft (x,y)/(r)d^(y) = f(X), v/ e w h , vx e n. 

Our second domain will be the subset Q aolm = {Z e C NxN : Z* Z = ZZ*} of 
all normal matrices in C NxN . By the Spectral Theorem, any Z 6 £! norm can be 
written in the form 

(2.5) Z = U*DU, 

with U € U(N) unitary and D diagonal; D is determined by Z uniquely up to 
permutation of the diagonal elements, and if the latter are all distinct and their 
order has been fixed in some way, then U is unique up to left multiplication by 
a diagonal matrix with unimodular elements. Consequently, there exists a unique 
measure dfih(Z) on f2 norm such that 

f(Z)d t i h (Z) = (irh)- N [ f f{U*DU)e-\\ D \\ 2 ' h dUdD V/, 

JU(N) JC N 



2 The authors are grateful to M. Bertola for this result. 
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where dU is the normalized Haar measure on U (N) , and dD is the Lebesgue measure 
on C", where we are identifying the diagonal matrix D — diag(di, . . . , d/v) with 

the vector d = (d i,...,d N ) S C", and ||D|| 2 = ||d|| 2 := |<ii| 2 H h \d N \ 2 . Again, 

one easily checks |AEG| that 

(2.6) f Z*iZ k dii h {Z) = 5 ]k k\h k I, 

so that the elements ()2.3|) are orthogonal also in L^, N (f2 norm , dun), and we let Tit 
be the subspace spanned by them. The reproducing kernel is then given by 

(2-7) tf fc( x,Y) =E 

fe=0 

(with the series converging for all A, Y G C NxN ), in the sense that 

K h (X, Y)f(Y) dfi h (Y) = f(X), V/ e H h , VA e Q noim . 

!l norm 

Remark. At first sight, the most natural candidate for the vector- valued space 
TL would seem to be the subspace L? ol CN (fl, dfi) of all holomorphic functions in 
L^, N (fl, d[i) (i.e. of all square-integrable C^-valued functions which depend holo- 
morphically on the coordinates z±i, . . . ,znn of the point z € f2). However, this 
choice turns out to be too simple-minded: the reproducing kernel is then just 
k(x,y)I, where k{x,y) is the reproducing kernel of the ordinary (scalar-valued) 
space L^ ol (fi, d/i); the Toeplitz operator (to be introduced in the next sec- 
tion) is just the N x N matrix [T,p jk }^ k=1 of Toeplitz operators on L^ ol (Q, dfx); 
and the Berezin transforms (also to be introduced in the next section) are just 

T4, = p0 Jfc ]f fe= i and T T V , = [EiLi T ^i T i>ik]f,k=v Thus > for instance, for the 
spaces (H7J with n = 1 (i.e. on fl = C), \\T { J l) T^ ] - T^\\ -> as h 0, while 

\\2ir\rp(h) rp(h)i _r r (h) || n 

II m v 1 ^ > 1 i> \ I4>,ipl" ' 

where 

\<t>, tjj] ■— d<j> ■ dip — di\> ■ d<j) 

(here d, d are applied individually to each element of a matrix, and the dot stands 
for matrix multiplication). From the physical point of view, this "matrix- valued 
Poisson bracket" seems to be a rather doubtful object, indicating that the spaces 
L? , CN are probably not the right route to take. 



3. The Berezin transform 

Let, quite generally, TL be a reproducing kernel subspace of C^-valued functions 
in L^-, N (fl, d/x), for some domain ft and measure [i on it, with reproducing kernel K; 
thus if is a C A ' XJV -valued function on ft and 

(3.1) f(X) = f K(X, Y)f(Y) dn(Y) VIefl, feH. 

Jn 

In particular, for any \ € C , the functions 



K Y , X {X) :=K(X,Y) X 
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belong to H, and 

(3.2) (/, K y , x )h = I K* Y , x f = X *f(Y). 

Jq 

See jMSl for more information on such spaces and their reproducing kernels. 

Let T be an arbitrary bounded linear operator on TL. For any fixed X G fi, 
the expression (TKx, x , Kx, x ') is evidently linear in x an d x'*\ thus there exists a 

unique N x N matrix T(X) such that 

X '*f(X)x=(TK x<x ,K XtXl ). 
We define the Berezin transform T of T by 

T(X) := if(X,X)- 1/2 f(X)if(X,X)- 1/2 . 

That is, 

x'*T(X)x = {TK XK ( XX) -i/-2 xl K XK{xx) -i,2 x ,) n - 
At first sight, this definition may seem a little ad hoc; the reason behind it is that 
this seems to be the only way to make the following statements true. 

Proposition 1. T is a C NxN -valued function on satisfying 

(i) T* = (?)*; 

(ii) i/ is a matrix-valued function on £1 suc/i i/iai <j)K x ^ x G 7i /or aZ/ X G 
and x G C , i/ien 

m^PO = #(x,x)- 1/2 (/>(x)ir(x,x) 1/2 

w/iere M^,/ := </>/; 

(iii) in particular, I{X) = I \/X G f2; 

(iv) ||f pOHc^c" < HTllw-^w, vx g a 

// iAe elements of Ji are holomorphic functions, then also 

(v) T(X) = VX only ifT = 0. 

Proof, (i) is immediate from the definition, while (ii) follows from the reproducing 
property l|3.1[) . and (iii) is a trivial special case of (ii). To prove (iv), observe that, 
by (22), 

\\Kx, x \\ 2 n = (K XiX ,K x , x ) n = x*#x,xP0 = x*OT,X) x - Ui^X,*) 1 /^^. 
Thus, for any x,x' G C", 

\x'*f{X)x\ = \(TK K{X j C) -i / 2 x ,K K(xx) -i / 2 x ,)\ 

< ll T ll II^X,_R-(A,A)-i/2 x || ||^(X,X)-i/2 x '|| 

= imi iixiic- k'iic«, 

and the assertion follows. _ 

Finally, for (v), note that T = implies T = 0, i.e. 

(3-3) (TK YtX ,,K XtX ) =0 V X ,x' 

whenever X = F. If 7Y consists of holomorphic functions, then K(X,Y) is holo- 
morphic in X and conjugate-holomorphic in thus the left-hand side of 1)3.3(1 is 
holomorphic in X and conjugate-holomorphic in Y. It is well known that such 
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functions arc uniquely determined by their restriction to the diagonal X = Y: 
consequently holds for all X, Y, i.e. 

X *(TK YtX ,)(X)=Q v x ,x'vx,r. 

Hence TKy, x > — for all Y and x' i an( i thus for any / € H 

x'*(T*f)(Y) = (T*f,K Y , x >) = (/, TKy, x ') = 0, 

i.e. T*f = 0. Thus T* = and T = 0. □ 

In analogy with the scalar- valued situation, we next define for any <j> £ C^ NxN (f2) 
the Toeplitz operator on 7i by the recipe 

T^f(X) := f K(X,Y)cf > (Y)f(Y)d f i(Y), 
Jn 

or, equivalently, 

T / = P(0/), 

where P : T 2 CN (f2, d/i) — > 7i is the orthogonal projection. Note that the last 
formula implies that HT^H^^-h < ||0||oo := sup^g^ ll < K-X')llc JV -+c Jv ; m particular, 
by Proposition^ also \\T<p(X)\\ c n^ c n < ||0||oo- 

Proposition 2. The following formulae hold: 

f+(X) ^K{X,X)- 1 ' 2 - [ K(X,Y)(j)(Y)K(Y,X)d^(Y) ■ K(X,X)-^ 2 ; 



T <p T^(X)=K{X,X)-^ 2 

K(X, Y)(j)(Y)K(Y, Z)ip(Z)K(Z, X) dfj,(Y) dfi(Z) 



>n Jn 

■K(X,X)-^ 2 . 

In particular, if (j) is a multiplier ofH {i.e. 4>f £H whenever f € TL), then 

f (X) = Kix.xy^^Kix.x) 1 ' 2 - 

and, similarly, when <p and tjj* are multipliers ofTL, then 

f^%(X) = AT 1/2 0WW,XMW(^^ 1/2 ■ 
4. Bad behaviour: all matrices 

We now exhibit an example of some pathological phenomena, showing in par- 
ticular that the straightforward generalizations of the expansions and (|1.12() 
cannot hold: first, apart from the integer powers of h, we will see that also \fh en- 
ters the picture; and second, instead of evaluations at X we get also contributions 
from other points. This applies to the case of the full matrix domain i\ = C NxN ; 
for the normal matrices f2 n orm, we will show in the next section that at least the 
second of these pathologies still remains. 

Theorem 3. Consider the full matrix domain SI ~ q NxN w ah N = 2. Let Tth 
be the spaces from Section\^ with reproducing kernels Kh given by \IH.4l - Let X be 
the matrix 

x= (o 0, 

Then: 
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, 01l(O)+022(O) n 

(i) ]hxiTr(X) = 2 u 
V 1 h^o * v ' \ 22 (O), 

Note that the matrix on the right-hand side does not depend on the value 

of <j) at X , but rather on its value at 0. 

(ii) If 4>(Y) = \ply22I , where y 22 denotes the (2,2)-entry ofY, then 



Tf- ] (X) = xVh + 0{h) ash^O. 



Proof. Since X 2 = 0, the series (|2.4|) for Kh(X, Y) becomes simply 
(4.1) 



XY* XY 
Kh(X, Y) = I -\ — = I 



c\h 



2h 



Thus 



T$ l) (X) = I K h (X,Y)4>(Y)K h (Y,X)dfi h (Y) 



(4.2) 



4(Y) 



XY*cj)(Y) (f>(Y)YX* XY*<j>(Y)YX* 



2/i 



2/v 



Ah 2 



dn h (Y). 



On the other hand, the change of variable Z = Y j\fh and Taylor's expansion imply 
that for any C°° function / on f2, 



(4.3) 



f(Y)d l x h (Y) = / f(VhZ)dm(Z) 



E-Uv(o)^, 



3=0 ■ 



where A = Ylj k=i d 2 /dzjkdzjk is the Laplacian on C 2x2 . Applying this to l|4.2J) . 
we therefore get 



X[A(Y^(Y))(0) + O(h)} 
2 

[A(4>(Y)Y)(0) + O(h)]X* 



X 



A(Y*<f>(Y)Y)(0) 



Ah 



0(1) 



X*. 



Now T^ l) (X) = K(X,X)- 1 I 2 T^ ) (X)K(X 1 X)- 1 / 2 - note that 



K h (X,X) 



-1/2 _ 



2h 



2/i+l 

1, 



10 S. TWAREQUE ALI, M. ENGLIS 

and thus, for any matrix A — ( Ql1 &12 

\a 2 i a 22 



H/2^/rrH/2-( 2h+i a ^ V 2h+i a ™ 



K(X,X)- L ' 2 AK(X,X) 



2h „ 
2iI+T a 2i a 22 



(4.4) 



o o 

1h 



K(X, X)- 1/2 AX*K(X, X)- 1 



-1/2 



2h + 



T a 12 0> 



21 1 



K(X, X)- 1 / 2 XAX*K(X, x)- 1 ' 2 = f 2^+1 a22 ° 



2h+T a ^ °j 

21 1 



Consequently, 



T * W -U(0)V2ft 22 (O) J + ° W 



o v|iA[r*0(r)] 22 (o) 





0(h) 



2MA{cf>(Y)Y} 22 (0) o) + ° {h) 



0(h). 



^A[Y*cf>(Y)Y} 22 (0) 


Note that, by a simple calculation, 

(4.5) A[r>(y)] afc (o) = £^(o), 

fe =i oyfea 

2 

A[^(y)y] o6 (o) = ^^(o), 
A[y>(y)r] a6 (0) = ^ Qb -Tr0(o). 

Letting h \ 0, (i) therefore follows immediately. For 0(F) = y 22 V%I, only the 
second term in the last formula for T^ h \x) gives a nonzero contribution, equal to 

lf\X) = (° Q ^^V22p22{0)^ + Q[h) = VJix + 0{hh 

as A[y 22 Y*] 22 (0) = 1 by g3J. This settles (ii). □ 

Remark. We pause to observe that applying (|4.3|l to f(Y) = Y* k Y k and comparing 
with 1|2.1[) shows that 

(4.6) A k (Y* k Y k )(Q) =k\c k I. 

This gives, conceivably, a way of evaluating the numbers without recourse to 
random matrix theory, and can also be used to show that the Ck have an interesting 
combinatorial meaning. Namely, expanding A fe and Y* k Y k yields 
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N 

A (Y Y ) = y ' ^iiji^hji ■ ■ ■ dikjk^ikjk 

h,ji,---,ik,jk = l 
N 

fc+iafc + 2 • • ■ Va 2 k-ibj 

ai,...,02is-i=l 

where, for the sake of brevity, we temporarily write dij for d 2 /dyij, and similarly 
for d^. Clearly a nonzero contribution only occurs if to each y there is applied 
precisely one d, and to each y precisely one d. (We also see that the result will be 
independent of Y, i.e. a constant.) Thus 

JV N 

/\ k (Y k Y k )= ^2 ^V(i)Jr(i)2/aia • ■ ' diT(k)3T(k)Va. k ak-i 

i 1 ,j 1 ,...,ik,jk = l ai,...,tt2*-i=l <y,re&k 

' ^*<r(l)i<r(l)J/«»fc<»fc+l • ■ ■ div{k)jv{k)ya.2k-lb- 

Changing the order of summations and using the fact that d^yu = Sik$jl, this 
becomes 

N N 

A (Y Y ) = 5^ 5^ 5^ ^ t (i)oi^jV(i)o ' • • • ^M*)a*^jV(*)»fc-i 

o-,reSfc ii j' fe = l ai,...,a 2 fc-i = l 

' ^<T(l)afe^<T(l)afc + l ' ' ■ ' ^Mfc) a 2fc-i < ^<7(fc) fc - 

Replacing i = . . . ,»*,), j = . . . ,j k ) by io t -1 J o r _1 and setting ^ = T _1 cr, 
we thus get 

' ^V(i)*fc^V(2)j M (i) ■ ■ ■ ^vM-vcfc-y^Vcfe) 6 ' 

In other words, [A fc (Y* fc Y fc )] a b//c! = Cfc5 a {, is the constant equal to the number of 
triples (//, i, j), where /U € 6fc and i, j G {1, . . . , iV} fc , such that 

j = (a,h, . . . ,ifc_i), j°M=(V(2)r'-.V()i)i il ). and *fe = V(i)- 

It is evident that, indeed, this number is zero for a ^ b (since the sequences 
(a, ii, . . . , ik~i,ik) and (wi), V(2) , ■ ■ ■ , V(fe)> ^) mus t be permutations of each other), 
while for a = 6 it is independent of a. It is also clear that Cfc is always an integer, 
a fact definitely not apparent from l|2.2[l . □ 

We conclude by giving a formula analogous to part (i) of the last theorem also 

for T^ h) T^ h) (X). It shows, in particular, that T^T^ (X)-T^ (X) need not tend 
to in general as h — > 0, but does so for scalar-valued <f> and ip. 

Theorem 4. Under the hypotheses of Theorem\^ 

/ Tr0(O) TrV-(O) n 
WtW/y^ 9 9 U 
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SI Jn 



Proof. Using the formula for T^T^ from Proposition [3 and 14.1(1 . we have 

T W T W (X)= f f ^ I+ ^^ { Y)K h (Y,ZMZ)(l+^-)d^(Y)d^ h (Z). 

Making again the change of variable Y i— > Yy/h, Z Z\/h, the double integral 
becomes 

v Z + ^^)^ Y ) K ^ ZMVhZ)(l + ^jf-) dfir(Y) d^(Z). 

The rest of the proof proceeds in the same way as in Theorem [21 using instead 
of H4.3JI the expansion from the following lemma (applied also to Y*<f>(Y) and ip(Z)Z 
in place of 4>(Y) and ip(Z), respectively), the fact that 

{{Y*<j>(Y),ip(Z)Z}} 22 = iTr0(O)TrV(O), 
and the relations l|4.4|) . We leave the details to the reader. □ 
Lemma 5. For (f), tp € C^ NXN (Q,), 

<j){v f hY)K l (Y, Z)ij){y/hZ) dpL X (Y) d^(Z) 

- 0(0)^(0) + h[Acf>(0) ■ yj(0) + 0(0) • AV(0) + {{0, ^}}(0)] + 0(h 2 ), 
as h — > ; where 

where E ik is the matrix [E ik ] ab = S la S kb . 

Proof. Using the Taylor expansions for <fi and ip, we see that the integral asymp- 
totically equals 

E lil+lffliliffl d a d P (j){0) f f g _s d~<d S ip(0) 
h 2 / y a y p K 1 {Y,Z)z 1 z d dm{Y)dm{Z) — i-A 

(The summation extends over all multiindices a,/?, 7,(5.) Note that the kernel 
satisfies Ki(Y,Z) = Ki(eY,eZ) for any e e C of modulus one; hence the last 
integral vanishes unless \a\ + \*y\ = \j3\ + \S\. Thus the coefficients at half-integer 
powers of h in fact vanish. The coefficient at h° is clearly 0(O)-0(O), since 



n Jn 



K 1 (Y 1 Z)dm(Y)d[ii(Z)= f Idm(Y)=I 

Jn 



>nJn Jn 
by the reproducing property of K\ and i|2.1|) . For the coefficient at h 1 , the only 
nonzero contributions, by virtue of the last observation, come from \a\ — \(3\ = 1, 
or I7I = \5\ = 1, or \a\ = \5\ = 1, or \J3\ = |7| = 1 (i.e. from yy, zz, yz, ovyz). Since 



yijVki K i( Y , z ) d(J-i(Y) d(ii(Z) = / yijy k il dm(Y) = SikSjiI 
inJn Jn 

(and similarly for ZijZki), the first two possibilities contribute 
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and 4>(0) ■ Aip(0), respectively. For the yz possibility, the corresponding integral 
vanishes, since the integrand is a function holomorphic in the entries of Y and Z* 
and vanishing at the origin. Finally, for the last possibility yz we use the series 
(12. 4|) to split the integral as 



(4.7) / / y ij z kl K l {Y,Z)d l JL 1 {Y)dni{Z) 

Jn Jn 

= E — ( / V tJ Y m dm(Y))( I z kl Z* m d^(Z) 

Using again the invariance of <i/ii under the change of variable Y i— ► eY, we see 
that we only get nonzero contribution for m = 1. In that case, 

VijYdfj,i(Y) = / y ij yabdni(Y) = 5 ia Sjb = [Eij]ab, 

and similarly for the Z integral. Thus the integral Q4.7H equals 

—EijEik = -SjiEik, 

and the total contribution from the yz possibility is 

dip 
dzu 



i,j,k,l 

which concludes the proof of the lemma. □ 
5. Bad behaviour: normal matrices 

The matrix X = ^ ^ featuring in the last section was not normal; this 

tempts one to hope that things might perhaps still work out fine for the domain 
fi norm of all normal matrices. We show that even in this case, unfortunately, the 
non-local behaviour described above still persists. 

Theorem 6. Consider the domain fi norm of all normal N x N matrices, with 
N = 2. Let TLh be the spaces from Section^ with reproducing kernels given 
by \2.7$ . Let X be the matrix 

of the projection onto the first coordinate. Then 

[Tf\x)] 11 = [Tf\l)] 11 , 

[Tf ) (X)] 22 = [T^(0)] 22 . 
Consequently, the asymptotic expansion (1 1.1 1\ cannot hold. 
Proof. As X is a projection, we have X J = X Vj > 1; thus 
x \ .V 'V •' x \ Y* j 

K h (x,Y)=j2-^= I + x zZii h i = ( I - x )+ XK ^n 

3=0 J ' j=l J ' 

Thus 

T^\x) = f K h (X,Y)<f>(Y)K h (Y,X)dti h (Y) 



14 S. TWAREQUE ALI, M. ENGLIS 

= X- [ K h (I,Y)cj)(Y)K h (Y,I)dn h {Y)-X 
K h {I,Y)cj>{Y)dnh{Y) -(I-X) 
4>{Y)K h {YJ)d^ h {Y)-X 



+ x 


/ 




Jn DOI 


+ (I 


-X). 


+ (I 


-X). 


( an 


a\ 2 \ 


1^21 


»22 / 



But for any matrix A = ( ^ ) , we have XAX = ( °jj x jj ) , XA(J - X) 



a 12 




, etc.; hence 

[Tf\X)}n = [jf ^(7,y)0(y)^(F ! 7)d/i,(y)] ii = [Tf\j)]n, 
and similarly, since Kh(Q,Y) = I, 

[Tf\x)h2 = \ I 4>(Y)dv h (Y)] = [Tf\0)] 22 . 

Finally, since K h (X,X) = {^^ , K h (I,I) = e x l h I, and 1^(0,0) = I, the 

assertion about rj ft) (X) = K[X, X)-^ 2 T^ h \x)K(X, X)- 1 ' 2 follows. □ 

For completeness, we also state the analog of Theorem 0J which shows, among 
others, that the expansion 1)1. 12[) cannot hold. Its proof is the same as for Theo- 
rem 

Theorem 7. In the situation of the preceding theorem, 

[Tf^\x)] u = [T^T^\l)] lu 
[T { J l hi h \x)} 22 = [T { J l) T^(0)} 22 . 

6. An application of stationary phase 

In this section we finally start exhibiting also results in the positive direction, 
namely, by using the stationary phase method we establish the existence of an 

(albeit non-local) semiclassical asymptotic expansion for Tr^ for the case of the 
normal matrices. 

Recall that the stationary phase (WJKB) method tells us that if S, 4> are smooth 
complex- valued functions on some domain in C" , such that S has a unique critical 
point xq (i.e. S'(xo) — 0), which is nondegenerate (i.e. detS"'(xo)) ^ 0) and is a 
global maximum for Re S, and 4> is compactly supported, then the integral 



(6.1) 



h~ n / (j>{x) e s ^/ h dx 
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has an asymptotic expansion 

oo 

(6.2) e s(,*o)/h h? £j0(x o ) as h -> 0, 

with some differential operators Cj whose coefficients are given by universal expres- 
sions in S and its partial derivatives. See e.g. | Hrm| . Section 7.7. The hypothesis 
of the compact support of cj> can be replaced by the requirement that the integral 
(|6.1(l exist for some h = ho > 0, and that the maximum of Re S at xq strictly 
dominate also the values of Re S at the boundary or at infinity, in the sense that 
Re S(x n ) — + Re S(xq) =P- x n — ► xq. 

On the other hand, if the global maximum of ReS 1 is not a critical point, then 
(|6.1(l decays faster than any power of h as h \ 0. 

The formulas for the operators Cj are fairly complicated in general, but fortu- 
nately become quite explicit if the phase function S is quadratic (which will be the 
only case we will need). Namely, assume that 

S(x) = -{Q(x - x Q ), x - a; ) C " 

for some matrix Q with positive real part. Then xq is a unique critical point of S, 
is nondegenerate, and 

(6.3) Cj = ^Q s , where Q--(Q- 1 9,9). 

Let us now apply this to the integral defining T^ h \x) in the case of the domain 
of normal matrices, viz. 

Ti h \x) = [ if /l (x,x)-V2^ /i (x,y)0(y) J Fr h (y,x)^(x,x)- 1 / 2 d Mh (y). 

Let 



(6.4) Y = UDU*, X = VCV* 

be the spectral decompositions of Y and X, respectively. Observe that owing to the 
invariance of the kernels Kh and the measures \ih under unitary transformations, 

we have T^ h) (X) = VT ( p{C)V*, where cjj v (Y) := V*(j>{VYV*)V; thus it suffices 
to deal with the case of V = I, i.e. when X = C = diag(ci, . . . , cat) is a diagonal 
matrix. From (|6.4|l . we then have 



[K h (X,Y)] ij = 



°° [C k UD* k U*] % 



^ klh ! 

k=0 

oo N k ~l k * 
\ - \ - c i Uil a l V-lj 

fe=0 1=1 
N 

J2e Cidl/h UiiUji, 
i=i 
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and similarly for Kh(X, X). Thus the matrix entries of Ti (X) are given by 
— r r N 

[TP{X)] ab = {nh)- N / V e-^ 2 / 2h e c ° d '/ h u al u jl( f> 3k (UDU*) 

JU(N) J C « j } k,l,m=l 

■ e^ h u bm u km e^ 2 / 2h e-^^ h dUdD. 

For simplicity, we will write (j>(U ; d±, . . . , dfj) instead of <fi(UDU*). The last integral 
over D is precisely of the form l|6.1[l , with phase function given by 

S(di, . . .,d N ) = c a di + c b d m - \\d\\ 2 - ^ ^ ^ . 

The critical point condition S' = amounts to 

Ca^ii = di, c b S ml — di, Vi = l,.,.,N. 

It follows that there is no critical point if c a ^ c b , or if c a = c b and I =/= m; 
while for c a = c b ^= and I = m, or c a = c b = and I, m arbitrary, there is a unique 
critical point 

l-th slot 

d = (0, . . . , 0, c a , 0, . . . , 0) = c a xi , 

which satisfies the assumptions for the application of the stationary phase method. 
The critical value is 

Cf \ I 1 2 i I 1 2 I 1 2 \ c a\ 2 + \c a \ 2 „ 
d{CaXl) = \Ca\ +\Ca\ ~ \C a \ = 0, 

and the operators Cj are equal to j*^, by l|6.H[l . By ({6.2(1 . it therefore follows that 
[Tf-\x)] ab = 0(h°°) for c a ^ c b , while 

N oo , r „ 

(6.5) [T^ ] (X)] ab ^ Uaiu^uuUu (Al d) <f> jk )(U;c a xi)dU 



j,k,l=l r=0 

as h ~ * if c a = c b ^ 0, and 

JV oo 



[T^ h \x)] ab w — / UaiUjiu km u bm (A r {d) (j) jk )(U-,0)dU 

j,k,l,m, = lr=0 r ' ^ U ( N ) 

= E 77 / (Afa0a6)(tf;O)dtf 



7 

(6.6) = ^ — (A[ d) Qb )(O) (as 0(C7; 0) = 0(0) is independent of £/) 



r 

r=0 



as /i — > if c a = c b = 0. Here the subscript at A indicates that it applies only to 
the d-variables in 4>(U; d\, . . . , d^). 

Thus the coefficients at each h r in the asymptotic expansion do not depend on 
the jet of (f> at X, but rather on the behaviour of near the whole orbit {UP a U* : 
U € U (N)} of the spectral components P a := diag(0, . . . , 0, c a , 0, . . . , 0) of X. Also, 
the off-diagonal entries asymptotically vanish (i.e. are 0(/i°°)) if c a ^ c b , which is 
quite unexpected. 
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Observe that setting c a — in Ij6.5(l gives 

- — - N oo , r 

E E h ^- A w^(°) 

where 



r! 



Kjk ■= E / UalUjlUklUbl dU. 

7=1 J U(N) 



It can be shown that 

„x S a jS kb + S ab 5 k j 
(6.7) - — 



3 N + 1 



and in fact, 



(6-8) / UaiUjiUkiUM dU = — . 

(The above relation can be obtained by an application of standard orthogonality 
relations for the matrix elements of irreducible representations of compact groups 
- in this case applied to the irreducible subrepresentation of U (N) , carried by 
second order symmetric tensors, in the decomposition of the natural representation 
of U(N) ®U(N).) Thus we have 



K ] u « E 



r=0 



W A r (d) (^ ab + 6 ab Tr (A)(0) 
rl N+l 



which is different from l|tj.6f) . Thus we see that, in general, it is not possible to use 
the formula l|6.5(l in both cases c a = c b ^ and c a = c b = 0. 
For scalar-valued <fi, Ijfi.SJl simplifies to 



oo N 



[Ti h \x)] ab « EE / ™ (A r {d) <i>)(U;c a xi)dU; 

r=0 1=1 ' 

and if in addition is independent of U, i.e. (p(UYU*) = 0(F) W S E/(iV), then 
the last integral can be evaluated by Schur's orthogonality relations, yielding 

- — - oo , r 

(6-9) [Tf\X)] ab « Sab ■ E 7T ( A ( d )^)( c «Xa). 



r=0 

In the general case, however, it does not seem that ()6.5|) can be simplified in any way. 
In the same manner, one can also prove the following formula for the asymptotics 

oiT^T^, which of course reduces to Ij6.5|) upon taking for ip the constant function 
equal to /. The strange-looking operators A4 mq originate from the formula 

Theorem 8. For any functions € C^? N y N (£l noTm ) and a diagonal matrix X — 
diag(ci, . . . ,c N ), 

[Ti h) T^\x)] ab = 0(h°°) ash^O 
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if C a ^ Q>.' 

N oo 

1 1' I I _ _ _ 

U ( i' n T t 'U'j im Uj rn Up rn WpqWf i q 



: , ,. , n r - JU N) JV 



(6.10) m,p,q=l 



U(N) JU(N) 

WlqW bq (M^^ijlpkOiU, W\ C a Xm, C a \q) dU dW 



as h — > if c a — ci, ^ 0, where 
(M r mq 0M(U,W;d,e) := 
and 



(A (d) + A (e) + & 3 (U;d)^ kl (W;e) 

i T !t )T i h} (X)}ab~ V V / / UaLUiLUjmUpm 

(D.llj p,q,L,M=l 

WpqWkqWlMWbM (M^qfiijlpkOiU, W \ 0, 0) d[/ rfM^ 



as /i — > if c a = Cb = 0. 

The formula can clearly be simplified upon carrying out the summations 

over L and M and performing the two integrations (which can be done since 
(j>(U;0) = 4>(0) and ip(W;0) — ?/>(0) are independent of U and W) via Schur's 
orthogonality relations; the result is 

, * N oo , r 

[Tl h) T$\x)] ah K E ^AC g (#) a6 (0,0). 

m : q—l r—0 

Similarly, as with l|t).5|l and l|t).6|l . using (|6.7|) it can be shown that for c a = the 
formula (|6.10() reduces to 

A " h r X^[(0 + /Tr^)(y; + /Tr^)] afc (O,O) 
r! (iV+1) 2 



[TfTf(X)] Qi) « ^ E 

m,q— 1 r—0 



which is different from . 1 1|1 . We refrain from going into these details because they 
are not needed anywhere in the sequel. 

We conclude this section by observing that (|1.6() also fails in general. 

Proposition 9. Let <f> he the C NxN -valued function on fl n0 im, N > 2, defined by 
4>{Z) = |dctZ| 2 /. Then 

limTf'pO^O VXei2 norm . 

Proof. Since (j)(VZV*) = V<j)(Z)V* W € U(N), it is enough to check the assertion 
for diagonal X, so let X = diag(ci, . . . , c^r). As 4>jk{U ; d) = 5jk\di . . . d^\ 2 , we have 



(A( d )(j)jk){U;d) = Sjk }Jdx...d m ...d 

rn 



IN? 
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3jk 2 ^ \dl . . . d mi . . . d m2 . . . d]y\ 2 , 



(A r ^ jk )(U;d) = 8 jk r\ 2J \di . . . d mi . . . d„ l2 . . . d mr . . . d 



N\ 2 , 



mi<m2<-<m r 



(A[ d) <p jk )(U;d)=0 for r > N. 

(Here the hat " indicates that the corresponding variable is omitted.) Thus by (|6.9I) 
and (ETUI 

[T^U = SabdCafh"- 1 + h N ) + 0(h°°), 

that is, 

lf\x) = h N - x X*X + h N + 0{h°°) 
as h — ► 0, and the assertion follows. □ 

Similarly, it can be shown that breaks down too: for instance, for (f>(Z) = 
\AetZ\ 2 e~ T ' c( - z ' z h, one has 

(where Halloo := sup^gf; ||0(^)||c w ^c iv )- This can be proved by observing that 
the operator is diagonal with respect to the basis l|2.3[) . with eigenvalues 

{k+l)h N 
(h + l) 2JV + fe 5 

and sup fc (fc + l)/(h + l) fc « l/(e/i) = \\<j>\\H N /h. We omit the details. 

We now turn to classes of observables <fi which are more manageable than the 
general case. 

7. Spectral and [/-invariant functions 

A function <j)(Z) of Z e f2 norm will be called spectral if it is a function of Z in 
the sense of the Spectral Theorem: that is, if there exists a function / : C — > C 
such that <j> = where 

(7.1) f*(Z) :=C/-diag J (/(d J ))-C/* if Z = U ■ diag^) • U*. 

Our first observation is that for spectral functions, all goes fine with the Berezin- 
Toeplitz quantization. 

Theorem 10. If <p = f# and ip = g# are two smooth spectral functions, then there 
exist unique spectral functions p r> r = 0,1,2, ... , such that 



,(h)rpih) _ 

in the sense of operator norms (i.e. as in In fact, 

Pr = C r (f,g)*, 

where 

(7.2) C r {f ) g) = ~d r f-tfg 

r\ 
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are the operators tyl.Kfy for n = 1 . 

Proof. Recall that the monomials z k , k — 0,1,2,..., are orthogonal in the Segal- 
Bargmann space (|1.7I) for n = 1: 

(^V) LL](c ^ h) = S kl k\h k . 

Comparing this with (|2.6(l . we see that the mapping 

(7.3) i : Z k Xj ^z k ® Xj 

J hol<. 



is a unitary isomorphism of our space Tih onto the tensor product L^ ol (C, d/i^) 



C^. Now if <f> = f# is a spectral function and %, rj E C N , then 
(T^^ fc X ,^ry) = ^Z k X ,Z l r,) L , cN(n _^ h) 

V *Z* [ ( f ) (Z)Z k xd^ h (Z) 

dD 



f [ V *UD* l <p(D)D k U* X dUe- Tl ( D * D '> , x . 

ic« Ju(n) \^ h r 



(7.4) / UXU* dU = J. 



However, for any matrix X 

f UXU* dU = - 

- U(N) N 

(Indeed, performing the change of variable U ^ U\U and using the invariance of the 
Haar measure, it transpires that the left-hand side commutes with any U\ € U(N). 
Thus it must be a multiple of the identity. Taking traces and using the cyclicity of 
the trace, (|7.4() follows.) Thus we can continue the above calculation with 



JV 

-\\df/h. 



= {x,v)±Y,f c J j d*f(d jJ , {nh)N 
1 N 



(why 

dD 



N ■ 

3=1 

= { X ,V) (Ti h) z k ,z l ) LLi(c ^ h) 

= {{Tf ] ® I){z k ®x),z l ® v}LlJc.d^ h )®c« ■ 

Consequently, under the isomorphism l, the operator on Tih corresponds to 

the operator ® I on L^ ol (C,d/x^) ® C", and the desired assertions follow 
immediately from the ordinary Berezin-Toeplitz quantization on C. □ 

We list one more corollary of the above isomorphism l; it will not be needed in 
the sequel, but should be contrasted with Proposition|5]at the end of Section lU and 
the example immediately thereafter. We omit the proof. 

Proposition 11. For any spectral function <f> = f# and x G C, 
(7.5) Tf\xl)=jf\x)-I 
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where the Tj on the right-hand side is the ordinary scalar-valued Berezin trans- 
form of the operator Tj on L^ ol (C, d/i^). In particular, 

limj|lf ) || oo = lim||rf|| = ||0|| oo . 



Remark. We pause to note that for the full matrix domain ft — C , the spaces 
Tih are not isomorphic to L^ o] (C,dv) for any rotation invariant measure v on C 
The reason is that the numbers c& in (|2.2|l . which take over the role of the A;!, are 
not the moment sequence of any measure on [0, oo) if N > 1. This can be seen by 
checking that 

k+l 



3 = 1 

where 

N-l 



( » + ij5niyn(*+J)-/ W**w 



di/jv(2) := - > ^77 — k J e N dz; 



7T 

J=0 



thus for Cfe to be a moment sequence (even of a measure which is not necessarily 
non-negative) it is necessary and sufficient that 

(7 - 6) fo+T^n<"-fl H 

be a moment sequence. However, the latter cannot be the case, since (|T.f>l) has only 
a finite number of nonzero terms. 

We restrict our attention exclusively to f2 n orm in the rest of this paper. □ 

Returning to the main line of discussion, we proceed to introduce another class 
of functions. 

A C Nx ^-valued function <f> on f2 n orm will be called U -invariant if 

(7.7) 4>{UZU*) = U <j>(Z) U* Vf7 G U(N) VZ e ft norm . 

Clearly, a spectral function is //-invariant, but not vice versa: an example is the 
function (j>{Z) = |detZ| 2 / from the end of Section [SI The relationship between 
spectral and [/-invariant functions is clarified in the next proposition. 

Proposition 12. A function <f> is U -invariant if and only if there exists a func- 
tion f(d\;d2, ■ ■ ■ ,d]Sf) from C x C w_1 into C, symmetric in the N — 1 variables 
c?2, ■ ■ ■ j djy, such that <j) — f#, where 

(7.8) f*{UDU*) := U ■ dia gi (/(d i; d u . . . , d h . . . , d N )) ■ U*. 

The function f is uniquely determined by <f>. 

Further, <f> is spectral if and only if f depends only on the first variable, i. e. if 
and only if f{d 1 ;d 2 , ■ ■ -,d N ) = /(di;0, ... ,0). 

Proof. For any complex numbers ex,...,ejv of modulus one, consider the matrix 
e = diag(ei, . . . , ejv). Then e £ U(N) and eDe* = D for any diagonal matrix D; 
thus by (|77|) 

<f>{D) = e<t>{D)e* Ve u . . . ,e N e T. 
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Consequently, 4>{D) is also a diagonal matrix. Define the functions /i, . . . , /at on 
C^by 

(7.9) f j (d 1 ;d 2 , . . . ,d N ) := (j)jj(D) where D = diag(di, d N ). 

For any permutation a of the set {1, 2, . . . , iV}, let F a denote the permutation 
matrix [F a ]jk — 5 CT (j),k- Then F a £ U(N) and 

F a DF* = diag(d CT (i), . . . ,d CT (jv)) if -D = diag(di, . . . , djv). 

Thus by fT7| again 

fa{j)(di;d2, . . . ,d N ) — (rfo-(i) ; rf CT ( 2 ) , . . . , d a ^ N - ) ). 

It follows that /j is symmetric with respect to the N — 1 variables di, . . . ,dj, . . . ,dx 
and = / # for / = /i- 

Conversely, it is easily seen that any function of the form (|7.8(l is [/-invariant, 
and f* = g* / = g. 

Finally, the assertion concerning spectral functions is immediate upon comparing 
G3 and £□). □ 

One consequence of the last proposition is that the mapping 

(7.10) f# — > (/ b ) # 
with f : C C defined by 

/ b W :=/(z;0,...,0) 

is a projection from [/-invariant functions onto spectral functions. (Here the first 
# in H7.10I) is the one for [/-invariant functions from l|7.8|l . while the second is the 
one for spectral functions from Ij7.1|l : however, there is no danger of confusion in 
this abuse of notation.) In terms of = <f>, the function f b can be expressed 
directly by 

f\z) = <t>u(zE n ), zeC, 
where En is the matrix of projection onto the first coordinate, i.e. [-Enjjfc = SijSik- 
The projections / i— > f b and (|7.10(l will play a crucial role in the next section. 

8. Quantization of [/-invariant functions 

We now proceed to establish our final result — a generalization of TheoremllUlto 
[/-invariant functions. The key ingredient is played by the following specializations 
of the asymptotic expansions from Section [B] 

Theorem 13. For any smooth U -invariant functions <j) = andi/j = g# on Q nolm , 

, — . oo 

(8.1) Tj'"^^^ 

r=0 



and 

. s oo 

(8.2) r< k) f ^^.(W)* 

r=0 

as h — > 0, where l r (f> and m r ((j), tp) are the functions on C defined by 

(8.3) W(z):=i(A'7)(z;0,...,0) 
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(that is, l r <j) — ^(A!"„()i) b ) and 



2:S 



.4) m r (<j>,ip)(z) 



If. 9 



2 \ r 



d=(z;0,...,Q) 
e=(z;0,...,0) 



Proof. In principle this could be gleaned from the formulas H6.5J) and (|6.6(l . but 
it is better to use directly the definitions: if X = VCV* with V £ U(N) and 
C = diag(ci, . . . , cat), then by Proposition [2] 

Ti h \x)= f x /l (x,x)- 1 / 2 ^(x,y)0(y)^ h (y,x)^(x,x)- 1 / 2 rf/i h (y) 



/ / vK h (c,cr'i' 2 v*f^ vckv 2? ku ' kudu*) 

Jc« Ju(N) ^ k\h k 



^ Z!^ dUe ^ 

C k V*UD* k 



1=0 



f f ye-cc-^C^V^ 

Jc» JU(N) ^ k]h 



^ DlU *^ C * 1 c -CC-/2hy* dUe -Tt(D*D)/h dD 



llh 1 (nh) N 

=^L Vc - cc '' 2 ^^ T ' {D ^ {D)D ' } 



. e -CC'/2hy* e -Tr{D'D)/h 



kj 

dD 



(nh) N 

N 



- Y I Ve- cc '/ h e^ c+d ^/ h ^ n (D) V* e~^l h 
N ~{Jc« ( 7rh y 

( 1 N f 

V ■ diag fc [jr22 f(dj;di, ...,dj,...,d 



V ■ diag 



v - di ' A &k( E ^(A r /)(c fc ;0,...,0)) -V* 

r=0 



r=0 

Here we have used, in turn, the formula i|2.7[l for K/ l (X, Y); the [7-invariance of <fi; 
the formula J731 for the integral over U(N); the fact that Tr(D* k (f>(D)D l ) = 
Sj=i djd l j<f)jj(D) : combined with the summation of the exponential series and the 
commutativity of C with C* ; the fact that <j> — f# ; the independence of the integral 
on j; the stationary phase expansion; and (|7.1(l and the definition of l r <p. 
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The proof of Ij8.2|) is similar: 3 
T W T W {X )= [ f K h (X,Xr^K h (X,Y)4>(Y)K h (Y,Z)^(Z) 

■ K h (Z, X)K h (X, X)- 1 ' 2 dfi h (Y) dn h (Z) 



Ve -cc> /2h E C^riJD^ m E D l U*WE* 1 
C N JC N JU(N) JU(N) k k'.h l 



llh 1 



<P( E ) E .,,1 e~ cc *l 2h V* dUdWd i x h {D)d i i h {E) 



1 



1 



c N Jc N 



y e -CC*/2h ^ 



c k c 



k.L. 



l\k\m\h k + l + m 



Ti(D* k (j)(D)D l ) 



TY{E* l ^{E)E m )d^ h {D)d^ h {E) 



N 



c N Jc N 



J2 Ve- cc '' h e^ c+d ^ +& i c '^ h cl )ii (D)^ jj (E)d i x h (D)dii h {E) 



= V ■ diag fc 



N'- 



E 



g-| c fc| 2 / h e {d,iC k +diej+ejC k )/h 



C N JC N 



-(IMI| 2 + IM| 2 )/'* 



MD)^(E) 



dD dE 



(irh) N (irh) N 



■ V* 



V ■ diag fe 



g-|c fc | 2 /'! e (<iic fc +d 1 e 1 +e 1 - fc )/h e -(||d|| 2 + ||e|| 2 )/h 



= V • diag fe ( Y, 

oo 



f(d 1 ;d 2 , . . . , d N )g(ei, e 2 , . . . , ejv) 
9 2 



dD dE 



{nh) N (irh) N 



f(d)g(e) 



d=(c fc ;0,...,0) 
e=(c*;0,...,0) 



Here the penultimate line comes from the formula (|6.3() . 



□ 



Corollary 14. If g is a smooth function on C and g# t/ie corresponding spectral 
function on ri n orm, then 

(8.5) r 9 M:^ Ar 9) # ash ^°- 

In particular, 

(8.6) hm (X) = VX 5 = 0. 

Proof. Combine (|8.3|> with the last part of Proposition □ 



3 Of course, 18.11 can also be obtained from 18.21 upon setting xjj = I; but it is more instructive 
to give a separate proof. 
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The following theorem is the main result of this section and indeed, of this paper. 

Theorem 15. For any smooth U -invariant functions (f>,ip on Cl noTm , there exist 
uniquely determined functions go, gi, . . . , on C such that 



(8.7) T^T^^^hTT^ aah->0. 

m=0 

Moreover, if <f> = and ijj = g# , then the functions g m are given by 

(8.8) g m = Gm (f 1 g) b 

for some bidifferential operators G m onC N (independent of f andg). In particular, 

Go(f,g) b =/V, and 

(8.9) i 

Gi(f,g) b -G 1 (gJ) b = —{f,g b }, 
Zir 

the Poisson bracket of f b and g b on C. 

Proof. The uniqueness is immediate from (|8.(j|l . The existence is, by virtue of (|8.2|) 
and 18.5fl . equivalent to 



DC: 



r— m.n— 



Comparing the expressions at like powers of h on both sides, this becomes 

r ^ A n g 

71 — 

which is solved by the recursive recipe 

r i 

(8.10) g r = m r {^ rj>) - V -r A" 3r _„. 

* — ' n! 

From l|8.4[l it is also clear that g m are of the form l|8.8|l with appropriate bidifferential 
operators G m . Finally, a short computation using the special instances r = 0, 1 
of 

gives (Em. □ 

Remark. Note that the quantities G m (</>, ■0) b do no£ depend only on f b and <7 b : 
the bidifferential operators G m involve derivatives also in other variables than d\ , e\ , 
and only after these are applied one takes the restriction to di = • • • = d^ = e-i — 
■•• = ejv = 0. It is therefore quite remarkable that G\(4>, tl>) b — G\(ilj,(j)) b depends 
only on f b and g b — the derivatives with respect to the other variables having 
cancelled out. □ 

We indicate another proof of the last theorem, based on the isomorphism l|7.3|l . 
(We gave the proof above first since the isomorphism (|7.3|l is probably something pe- 
culiar to the domain of normal matrices, while the stationary phase method should 
work also in other situations. The proof below also requires a slightly stronger 
hypothesis on the functions 4> and ij>.) 
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For a function / on C N and h > 0, let Phf be the function on C defined by 

PkfM = I /(*, • • • , *n) c -(l*l"+-+l^l")/fc 

Theorem 16. Let = /*, ip = g# be smooth U -invariant functions on f2 n orm 
such that the partial derivatives of f and g of all orders are bounded, and let C r be 
the bidifferential operators Then 



T (h) T (h) _ > r rp(h) 

L 4> 1 i> ~ " 1 c r (p h f,p h 



r=0 



g) 



# 



in the sense of operator norms. Consequently, \8. 7\ holds for 
■9>™= E - ! ^y^(A'V) b -9 r (A" £ 5 ) b , 

j,fc,r>0, ■'' ' T ' 
j+k+r—m 

where A' denotes the Laplacian with respect to the last N — 1 variables z%,...,zn. 
Proof. By a computation similar to the one in the proof of Theorem I1UI for any 

(T<pZ" X ,Z l v)= I ifZ*^(Z)Z k X dn h (Z) 
Jfl„ m 

= / / V *UD* l U*^(UDU*)UD k U* X dU e~ Tl{D * D)/h 

Jc» Ju(N) frfr)* 

= / / v *UD* l <j>(D)D k U*xdU e- T *( D * D V h 

Jc*> JU(N) (™) 

(by the J7-invariance of <p) 
dD 



N J c n (7T 



(nh) 



N 



(XiV} IZrit, f d l J d k j f(d. ] ;d l ,...,d. ] ,...,d N )e- 



\\dflh _ dD 



N 



\ \ - 'I: I d\d k f(d i; d 2: ...,...,d N )e-^ ' h s 



■4 ,t- »/ , , , n _iiwii 2 /i, 

!c N 

= (X,»7) {z k Phf,z l ) L 2 [Cjlllh) 
= (X,V) (T^z\z l ) Llai{C4lXh) 

= <( T pS ® ® *)> Z ' ® ^>LLi(C^)®C«- 

Consequently, under the isomorphism (|7.3|) , the operator on W ^ corresponds to 

the operator T p h \®I on L^ ol (C, d[Af l )®C N ■ Thus by the ordinary Berezin-Toeplitz 
quantization on C, 

rp(h)rp(h) ^ rp{h) rp{h) j 

1 <t > 1 4, ~ 1 P h f 1 P h g® 1 

OQ 

~^h r T Cr[PhfPhg) ®I 

r=0 
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V h r T {h) 

Z^' L 1 C r (P h f,P h g) 



# 



(the last isomorphism is the one from the proof of Theorem 110(1 . This proves the 
first claim. The second part of the theorem follows upon inserting the expansion 

oo I j 
3=0 J ' 

which follows from the Taylor formula (or stationary phase), and taking Berezin 
transforms on both sides. (The hypothesis of boundedness of the derivatives of / 
and g is needed in order that the resulting expansion for C r {Phf,Ph9) converge 
uniformly on C N , and thus imply the convergence of the corresponding expansion 

for T Cr{P h f,p h g)# b y tne inequality ||T^|| < ||0||oo-) □ 



For two U- invariant functions <j> — ip — g&, define their "star product" (/>*ijj 
as the formal power series 

oo 
r=Q 

As usual, this product can be extended by C[[h]] -linearity to all 4>, ip G the 
ring of all power series in h with coefficients in the algebra U of all [/-invariant 
functions on fi norm . Alternatively, upon identifying <f> = f# € IA with /, we may 
view this as the star product 

oo 

/*<?:=5> r G r (/,<?)> 

on the algebra S of all functions f(di; d%, . . . , djy) on CxC^ -1 symmetric in the last 
N — 1 variables, which again can be extended by C[[/i]] -linearity to all /,j g <5[[/i]], 
the ring of formal power series with coefficients in S. If we extend to <S[[/i]] by 
C[[/i]]-linearity also the operators G r , then the extended star-product will still 
satisfy the relations l|8.9l) . Further, * is clearly associative, since the multiplication 
of operators is associative — both (0 * ip) * rj and cf> *(?/;* 77) originate from the 
asymptotic expansion as h — > of [T^T^T^]~. (However, in contrast to a 
genuine star-product, the function constant one is not the unit element for *.) 

The appearance of / and g b , and not / and g, in l|8.9|l means that the C* -1 
part of / disappears in the semiclassical limit h — * 0, and only the projection / , 
which lives on C, survives. As mentioned before, we are dealing here with a quan- 
tum system which has N internal degrees of freedom. This is made clear by the 
isomorphism l|7.3|) . since the tensor product space L^ ol (C, d^h) ® C N is exactly the 
Hilbert space of a single quantum particle, moving on the phase space C and having 
N internal degrees of freedom. The full set of quantum observables of this system 
include those which do not have classical counterparts. The interesting fact that 
emerges from our analysis is that, it is exactly those observables which are Berezin 
quantized versions of [/-invariant functions, that have classical counterparts. Since 
the internal degrees of freedom are purely quantum in this case, they do not survive 
in the semi-classical limit. 
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